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FOREWORD 


Volume  I of  the  Final  Report  deals  with  Ocean  Micro - 
.-tructure  and  Volume  II  deals  with  Laser  Sensor  Technology- 
Assessment. 


ABSTRACT 


The  horizontal  and  vertical  length  scales  of  ocean  fine  structure 
appear  to  correlate  to  the  corresponding  wavelengths  of  inertial  waves 
(internal  waves  whose  motion  is  dominated  by  Coriolis  forces).  A prion, 
there  are  four  possible  non-linear  mechanisms  that  could  form  ocean 
fine  structure  from  inertial  waves;  dynamic  breaking  (shear  instabilities), 
static  breaking  (inverted  density  profile),  transverse  steepening  (reson- 
ance) and  longitudinal  steepening  (self  steepening  of  waves  riding  on  their 
own  velocity  fields).  Each  mechanism  has  been  assessed  on  the  basis  o 
the  available  data.  Of  the  four  mechanisms  proposed,  dynamic  wave 
breaking  is  shown  to  be  the  most  probable  mechanism  by  which  inertial 
waves  form  fine  structure.  This  mechanism  is  distinct  from  that  of 
"billow  turbulence"  due  to  an  unstable,  horizontally  propagating  gravity 
wave.  An  inertial  wave  is  a vertically  propagating  transverse  wave. 

The  transverse  velocity  field  of  the  wave  is  a horizontal  shear  layer  and 
it  is  the  instability  of  the  shear  layer  that  is  the  proposed  source  of  fine 

structure. 


jEfiECEDIlt}  PAGE  BLAMUNOT  FIIMSD 


table  of  contents 


abstract 

I.  introduction 

II.  INTERNAL  WAVE  SOLUTIONS 

III.  COMPARISON  OF  FINE  STRUCTURE  DATA 
WITH  INTERNAL  WAVES 

IV.  CONCLUSIONS 
REFERENCES 
APPENDIX  A 


-V- 


LIST  OF  ILLUSTRATIONS 


i Fig.  1 Idealized  Ocean  Fine  Structure  3 

Fig.  2 Temperature  Gradient  Record  Reproduced  from  Gregg  Gox 
and  Hacker.^  The  plot  begins  at  215m  in  the  upper  left  hand 
' corner  and  ends  at  400  m at  the  lower  right.  The  total 

length  of  the  arrow  represents  10  m in  depth,  Z positive 


upward.  4 

Fig.  3 Coordinate  System  on  Rotating  Earth  8 

Fig.  4 Density  Profile  Corresponding  to  a Statically  Stable 

Internal  Wave  13 

Fig.  5 Density  Profile  Corresponding  to  a Statically  Unstable 

Internal  Wave  14 

Fig.  6 Comparison  of  Internal  Waves  with  Fine  Structure  16 

Fig.  7 Wave  Amplitude  Required  for  Instability  19 

Fig.  8 Inertial  Wave  Stability  22 

Fig.  9 Estimates  of  Self  Steepening  Times  25 

Fig.  10a  Woods  Postulate  of  Billow  Turbulence  29 

Fig.  10b  Proposed  Source  of  Ocean  Fine  Structure  29 

I 

Fig.  A-1  Internal  Wave  Regimes  A-10 


I.  INTRODUCTION 


1-8 


Field  observations  have  revealed  that  the  ocean's  density  profile 
possesses  an  internal  structure  consisting  of  a series  of  thin,  laminar  "sheets' 
of  high  density  gradient  separated  by  "steps"  of  only  moderate  gradient.  An 


idealized  "fine  structure"  is  represented  schematically  in  Fig.  1.  Only  Woods' 


1 2-8 
observation  reveals  the  more  distinct  step-like  structure.  Later  observations 


illustrate  more  irregular  profiles  v'ith  a fine  structure  of  height  H and  a super- 
imposed microstructure  on  a length  scale  small  with  respect  to  H.  One  of  the 


more  irregular  profiles  is  illustrated  in  Fig.  2 which  is  reproduced  from  Gregg, 
Cox  and  Hacker.  ^ The  recording  demonstrates  temperature  gradient  as  a 


function  of  depth.  The  "spikes"  in  the  temperature  gradient  correspond  to  the 
"sheets"  of  high  gradient  and  the  small  scale  structure  in  the  temperature 
gradient  represents  the  microstructure  which  is  actually  present  in  the  "steps" 
of  the  idealized  fine  structure. 

The  "step  height"  H and  the  lateral  extent  of  the  observed  structure 
varies  widely  over  the  range  of  mean  density  gradient  indicated  in  Table  I.  The 
observation  time  refers  to  the  time  spent  taking  soundings  at  various  lateral 
locations  whereas  the  point  persistence  time  refers  to  the  time  duration  between 
soundings  at  a given  lateral  location.  The  term  persistence  need  not  necessarily 
imply  that  the  same  structure  remained  for  the  duration  of  the  persistence  time. 
It  is  always  possible  that  a second  similiar  structure  occupied  the  same  point 
at  a later  time. 


Various  mechanisms  have  been  proposed  to  explain  the  observed 


9 10 

phenomenon.  These  mechanisms  consist  of  double  diffusion,  interleaving, 


11 


dynamic  breaking  of  internal  gravity  waves  due  to  shear  instabilities  and  the 


12 


static  breaking  of  internal  waves  due  to  the  associated  density  inversions. 
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This  paper  is  not  intended  to  be  a review  or  assessment  of  the  proposed 
mechanisms.  What  v'e  wish  to  do  is  to  investigate  the  possibility  that 
internal  waves  are  responsible  for  the  generation  of  fine  structure.  Due 
the  the  sparse  amount  of  data  available,  we  can  only  attempt  to  determine  a 
mechanism  which  is  at  least  not  inconsistent  with  data.  Hereafter,  the 
terms  consistent  and  apparent  correlation  shall  be  used  to  infer  that  the 

data  does  not  demonstrate  otherwise. 

The  internal  wave  solutions  are  briefly  reviewed  in  Sec.  II.  To  avoid 
ambiguities,  let  us  adopt  the  notation  that  inertial  waves  are  internal  waves 
in  which  the  Coriolis  forces  play  a dominant  role.  Inertial  waves  propagate 
"almost"  vertically  with  a frequency  related  to  the  earth's  rotational  fre- 
quency. On  the  other  hand,  gravity  waves  are  internal  waves  which  are 
dominated  by  buoyancy  while  Coriolis  forces  play  a secondary  role.  In 
Sec.  Ill,  the  length  scales  of  internal  waves  are  compared  to  those  of  the 
fine  structure  and  it  appears  that  only  inertial  waves  can  be  responsible  for 
fine  structure.  This  suggests  that  a mechanism  by  which  inertial  waves 
could  form  fine  structure  is  inherent  in  the  difference  between  inertial 
waves  and  gravity  waves.  One  such  mechanism  is  shown  to  be  the  shear 
instability  of  the  wave.  Inertial  waves  are  much  less  stable  to  the  Kelvin 
Helmholtz  instability  than  are  gravity  waves . Other  mechanisms  are  sug- 
gested which  could  transform  the  inertial  wave  into  fine  structure  but  they 
are  shown  to  be  considerably  less  probable  than  dynamic  breaking.  This 

strongly  indicates  that  shear  instabilities  are  the  mechanism  by  which 

11 

inertial  waves  form  fine  structure.  Previous  investigators  have  proposed 
that  shear  instabilities  lead  to  ocean  fine  structure.  However,  their  concept 
of  "billow  turbulence"  due  to  a horizontally  propagating  gravity  wave  is  quite 
different  from  the  present  concept  of  an  unstable,  vertically  propagating  inerti^ 
wave.  A physical  description  of  the  proposed  model  and  its  distinction  from 
"billow  turbulence"  is  given  in  the  Conclusions. 
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II.  INTERNAL  WAVE  SOLUTIONS 


The  mathematical  description  of  internal  waves  is  a classic 
oceanographic  problem.  For  wave  frequencies  near  the  earth's  rotational 

(13)  . 

frequency,  the  Coriolis  forces  must  be  included.  Phillips  gives  an 
account  of  the  wave  solutions  in  the  presence  of  Coriolis  forces  but  that 
solution  is  carried  only  as  far  as  obtaining  the  wave  dispersion  equation. 

The  wave  solutions  for  the  velocity  and  density  are  necessary  in  order  to 
relate  ocean  fine  structure  to  internal  waves.  Hence,  we  must  expand  on 
Phillip's  results.  The  details  are  given  in  Appendix  A and  only  the  highlights 

of  the  model  and  solutions  are  discussed  here. 

We  adopt  a locally  orthogonal  system  on  the  surface  of  a rotating 

earth.  The  earth's  rotational  frequency  is  Q/2  and  the  coordinate  system 
is  illustrated  in  Fig.  3.  The  momentum  equations  include  buoyancy,  vis- 
cosity, and  the  Coriolis  forces.  The  flow  is  assumed  to  be  incompressible, 
that  is,  the  convective  derivative  of  density  is  zero.  The  incompressible 
continuity  equation  is  used  to  close  the  set  of  equations  in  pressure, 
density  and  velocity.  To  obtain  a solution  to  the  dynamic  equations,  we 
assume  that  the  only  velocities  are  those  due  to  the  internal  waves  and  we 
shall  expand  the  density  and  pressure  into  a static  and  wave  contribution: 

/-iv 

0 = p + p (z)  + p (x,  y,  z,  t),  (i) 

os  w 

and 

p = p - p g ^ + P (z)  + p (x,  y,  z,  t)  , (c) 

o o s w 

where  the  superscript  ~ denotes  a dimensional  quantity. 
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Fig.  3 Coordinate  System  on  Rotating  Earth 


The  static  density  gradient  must  be  specified  "a  priori".  Hence 
we  assert  that 


S P 
s 


(3) 


where  p is  the  reference  density  at  z - 0.  The  value  of  3 is  expressed  in 
o _4  , 

units  of  inverse  length  and  has  values  between  10  /m  in  the  summer 

thermocline  and  10  ^/m  in  the  deep  sea.  The  density  gradient  is  related 

to  the  Brunt  Vaisala  frequency,  N,  by 


N 


(4) 


where  g is  the  acceleration  of  gravity. 

In  developing  internal  wave  solutions,  we  introduce  the  following 

dimensional  variables:  a,  the  wave  frequency;  k^,  the  horizontal  wave 

number;  k , the  vertical  wave  number  and  u , the  fluid  velocity  in  the 
V P 

horizontal  direction.  These  dimensional  quantities  are  used  to  define  four 
dimensionless  parameters: 


e = u k /a  (Wave  Amplitude), 

P H 

= 'g/k  (Static  Ap/p  Across  Wave), 

V 


5 :=  k /k  (Vertical/Horizontal  Length  Scale), 

V H V 


and 

Re  = "p  a/pk  ^ (Wave  Reynolds  Number). 

wo  V 


Under  the  assumptions  that  N is  a constant,  £ « 1 (linearization),  A « 1 
(Boussinesq  approximation)  and  » 1 (Inviscid),  the  dispersion 

equation  becomes 
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(5) 


where  [i  is  a unit 

A 

|U  becomes 


vector  -11  ihe  direction  of  wave  propagation.  For  6^  « 1, 


(i  = k + 0(6  i)-t-0  (6  j). 

' i r V 


Clearly,  for  6 « 1,  the  internal  waves  propagate  almo  st  vertically. 

The  connection  between  ocean  fine  structure  and  this  type  of  internal  wave 
is  yet  to  be  established  but  the  discussion  will  be  restricted  to  this  limit. 
The  simplified  dispersion  equation  becomes,  for  6^  « 1, 


~2  ~2  .2  2 

a = N 6 + n,  , 

V J 


(6) 


= 2n  radians/12  hours  and  u is  the  latitude. 
The  wave  solution  for  the  density  is 


where  - 0 since  , 0 


= (p  e “H/lc  ) exp  (i  1), 


(7) 


w o 


ind  the  wave  function  is  given  by 


~ ~~  , ~ , T I V 

f = a t + k X + k y -t  k 

^ X y ' 


where  k^  and  k are  the  wave  numbers  in  the  x and  y directions  respectively. 

In  ordir  to  obtain  solutions  for  the  fluid  velocity  m the  x and  y 
directions,  we  utilize  the  limits  suggested  by  Eq.  (6).  The  frequency 


ists  of  a contribution  due  to  gravity,  6^^,  and  an  inertial  contribu- 
tion, Q The  dominance  of  the  buoyant  term  gives  rise  to  gravity  waves 

whereat  the  dominance  of  the  Coriolis  forces  leads  to  inertial  waves. 


a cons 
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For  Gravity  Waves  (N  6^  >Q^),  the  solutions  become 


a = N6  . (8) 

V 

u = -i  (k  /k  ) u exp  (i  $),  (9) 

w X H p 


and 

V = -i  (k  /k  ) u exp  (i  ^ ), 
w y H p 


(10) 


where 


For  Inertial  Waves  ^ obtain 

a = , 

w p y H X rl 


and 


(11) 

(12) 


(13) 


An  arbitrary  sum  of  the  plane  wave  solutions  over  all  positive  and 

negative  wave  numbers  is  a solution  to  the  equations  of  motion  if  and  only  if 

e is  much  less  than  unity.  However,  due  to  the  transverse  nature  ofthe 

waves,  situations  do  exist  where  the  nonlinear  terms  are  identically  zero 

and  the  solutions  are  valid  for  arbitrary  e.  Two  such  solutions  correspond 

(14) 

to  the  single  plane  wave  and  the  Resonant  Triad.  Using  one  plane  wave 

and  choosing  the  time  and  horizontal  location  such  that  the  wave  amplitude 
is  a maximum,  we  express  the  total  density  profile  for  arbitrary  S, 
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(14) 


where  is  2TT/k^.  Eq.  (14)  is  illustrated  in  Fig.  4 for  £ = 1.  The  density- 
profile  illustrates  an  "order  one"  wave  superimposed  on  the  mean  density. 
However,  some  distortion  (steepening  or  breaking)  of  that  wave  must  occur 
before  it  can  explain  the  field  observations.  The  density  profile  for  £ = 2 
is  demonstrated  in  Fig.  5.  Regions  of  the  density  profile  are  statically 
unstable.  Such  regions  could  overturn  or  flatten  to  form  steps  separated 
by  high  gradient  sheets.  ^ ^ This  does,  however,  require  waves  of  very  large 

amplitude.  We  shall  seek  a mechanism  that  could  transform  waves  of 
much  smaller  amplitude  into  fine  structure.  It  will  be  shown  that  the  most 
probable  mechanism  is  the  dynamic  breaking  of  inertial  waves. 
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Fig.  5 Density  Profile  Corresponding  to  a Statically  Unstable 
Internal  Wave 
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III.  COMPARISON  OF  FINE  STRUCTURE  DATA 
WITH  INTERNAL  WAVES 


While  we  are  not  yet  able  to  demonstrate  that  internal  waves 
actually  generate  the  ocean  fine  structure,  we  have  several  constraints 
on  the  waves  which  must  be  consistent  with  observation.  Otherwise, 
they  will  never  be  capable  of  explaining  the  origin  of  the  fine  structure. 

We  will  now  compare  the  fine  structure  length  scales  with  those  of  internal 
waves,  and  then  isolate  the  most  probable  mechanism  responsible  for 
generating  the  fine  structure  from  the  internal  waves. 

If  internal  waves  cause  fine  structure.  Figs.  4 and  5 suggest  that 
X will  correspond  to  the  step  H.  Furthermore,  we  assert  that  X , the 

V ^ _ 

horizontal  wavelength,  will  correspond  to  the  lateral  extent  of  the  fine 
structure.  Therefore, 

5 ^ Step  Height  ^ _ 

V ~ Lateral  Extent 

The  field  observations  of  5 are  compared  to  the  internal  wave  solutions 

V 

in  Fig.  6.  While  there  is  insufficient  data  to  suggest  a definite  correlation, 
the  data  that  is  available  does  appear  to  correlate  to  inertial  w’aves.  This 
suggests  that  the  mechanism  responsible  for  generating  fine  structure  may 
be  inherent  in  the  difference  between  inertial  and  gravity  waves.  That 
mechanism  will  be  shown  to  be  the  dynamic  stability  of  the  wave. 

There  is  a multitude  of  mechanisms  that  one  could  consider  as  a 
candidate  for  transforming  inertial  waves  into  fine  structure.  Basically, 
they  can  be  classified  into  two  major  categories:  stable  and  unstable. 
Internal  waves  may  be  driven  unstable  by  a density  inversion,  more 
commonly  referred  to  as  a static  instability,  or  by  a dynamic  instability. 

A primary  example  of  a dynamic  instability  is  the  Kelvin-Helmholtz,  or 
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Fig.  6 Comparison  of  Internal  Waves  with  Fine  Structure 
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shear  instability.  Stable  waves  may  be  capable  of  generating  fine  struc- 
ture through  their  nonlinear  interaction  with  the  velocity  fields  induced 
by  the  waves.  This  mechanism  can  be  divided  into  two  classes:  transverse 
and  self  steepening.  Transverse  steepening  simply  refers  to  situations 
where  resonance  occurs:  two  or  more  plane  waves  interact  through  their 
transverse  velocity  fields.  Self  steepening  refers  to  the  influence  of  a 
higher  order  longitudinal  velocity  field.  Such  a velocity  field  would 
transport  fluid  in  the  direction  of  wave  propagation  and  cause  self 
attenuation  of  the  wave.  These  four  mechanisms  may  be  summarized 
as  follows: 

1)  The  dynamic  breaking  of  a single  wave  due  to  shear  instabilities. 

2)  Breaking  of  a single  wave  by  static  instabilities. 

3)  Transverse  steepening  (resonance)  by  several  non-parallel 
plane  waves. 

4)  Self  steepening  of  a single  plane  wave  by  a higher  order 
longitudinal  velocity. 

In  addition  to  the  mechanisms  listed  above,  one  could  consider  the 
simultaneous  effect  of  the  sum  of  any  two,  three  or  four  mechanisms. 
However,  this  may  be  considered  over-ambitious  at  this  point  in  the 
development  of  a fine  structure  model. 

Let  us  now  assess  each  of  the  four  possible  mechanisms; 

The  dynamic  stability  of  the  wave  is  measured  by  the  Richardson 
number  J, 


j - AE  / \ - b ^ (1  - e cos  z)/e^  cp^,  (16) 

^ / V ^ / V 

0 6 z / \d  z / 

o 

where 

cp  = a/N. 
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For  J > 1,  a shear  instability  is  always  stabilized  by  the  stratification 
whe-reas  for  0 < J < 1,  the  wave  is  unstable,  provided  that  the  Reynolds  num- 
ber is  sufficiently  large.  A negative  Richardson  number  corresponds  to  a 
statically  unstable  wave.  This  occurs  when  e > 1 and  was  illustrated  in 

Fig.  5.  Subsequently,  the  discussion  of  dynamic  stability  will  be  restricted 
to  e < 1. 

For  gravity  waves,  <P  = 6^.  If  e < l it  follows  that  J > 1 and  the 
waves  are  dynamically  stable.  We  will  now  show  that  inertial  waves  are 
dynamically  unstable.  This  basic  difference  in  dynamic  stability  will 
explain  why  fine  structure  apparently  correlates  only  to  inertial  waves. 

Inertial  waves  correspond  to  6^  < cp.  Hence,  the  Richardson 
number  is  less  than  unity  for  e < i.  The  critical  value  of  the  wave 
amplitude  necessary  to  induce  unstable  wave  motion  is  obtained  from 
Eq.  (16)  for  J = 1. 


e 

c rit 


= 6 N/n 

V 3 


(17) 


The  critical  values  of  the  inertial  wave  amplitude  are  illustrated  in 

Fig.  7.  Wave  amplitudes  between  5 x lO'^  and  0.  5 are  required  to  explain 

available  data.  Hence,  dynamic  instabilities  would  occur  at  lower  wave 

amplitudes  than  would  static  instabilities  (e«*l).  The  mechanism  would  also 

explain  why  fine  structure  seems  to  accompany  only  inertial  waves. 

The  Richardson  criterion  for  wave  instability  has  been  satisfied. 

However,  a Reynolds  number  constraint  must  also  be  satisfied  before  we 

can  state  that  the  inertial  waves  which  correspond  to  fine  structure  are 

unstable.  Since  an  inertial  wave  is  a vertically  propagating  transverse 

wave,  it  creates  a horizontal  shear  flow  of  velocity  u and  thickness 
- 1 P 

(ky)  . The  stability  of  such  a flow  is  measured  by  the  flow  Reynolds 

number  Re^,  where 


NKyswwwwwww'.iK . 
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Fig.  7 Wave  Amplitude  Required  for  Instability 
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(18) 


Re 


f 


D U 

o P 

« 

k 


The  transition  value  of  Re^  is  a function  of  the  shape  of  the  fluid  velocity 
profile.  For  pipe  flow  and  boundary  layers,  the  transition  Reynolds 

, 16  U 1 u 

number  is  2300  and  2700,  respectively.  Free  shear  layers  have  an 
inflection  point  in  their  velocity  profile  and  are  less  stable  than  pipe  and 
boundary  layer  flows.  The  transition  Reynolds  number  for  a free  shear 
layer  is  of  the  order  of  10  to  10  . Using  the  mor e re strictive  value,  we 
assert  that  the  flow  is  unstable  when 


Re^  > 10^, 

or,  using  the  definition  of  £, 


X 

V 


^ 100  (2TT)^  V 

e a 


(19) 


6 2 ^ 

where  V = 2 x lO”  m /sec  and  o = Qy  Instability  will  occur  at  a value  of  e 

corresponding  to  £ . . Using  £ . for  £,  Eq.  (19)  becomes 

crit  crit 


X 


V 


> 


20TT 
~~  1/4 

(g3) 


(20) 


In  addition  to  the  constraint  imposed  by  Eq.  (20),  the  entire  analysis 
requires  that  the  mean  motion  be  undamped.  This  was  expressed  as 


-20- 


II  - I rTCfrn " - r.;--. L,.., 


or 


Re  > 1, 

V.' 


The  constraints  imposed  by  Eqs.  (20)  and  (21)  are  illustrated  in  Fig.  8. 

All  fine  structure  data  corresponds  to  dynamically  unstable  waves. 

This  strongly  indicates  that  shear  instabilities  are  the  mechanism  by 
which  inertial  waves  form  fine  structure.  A mechanism  utilizing  shear 
instabilities  was  originally  proposed  by  Woods  and  Wiley.  However,  the 
present  concept  of  unstable  horizontal  shear  layers  driven  by  vertically 
propagating  inertial  waves  is  quite  different  from  Wood's  concept  of 
"Billow  Turbulence"  driven  by  gravity  waves  which  propagate  horizontally 
along  density  "sheets". 

The  second  mechanism  by  which  inertial  waves  may  form  fine 
structure  is  the  static  breaking  of  the  wave  due  to  a density  inversion. 

This  process  was  illustrated  in  Fig.  5.  Large  amplitude  waves  raise 
the  heavier  fluid  to  a position  above  the  relatively  lighter  fluid.  Static 
instabilities  then  break  the  wave  and  a staircase  density  profile  results. 
While  this  is  a reasonable  postulate,  it  fails  to  explain  why  fine  structure 
seems  to  correlate  only  to  inertial  waves.  In  addition,  and  most  funda- 
mental, the  mechanism  requires  wave  amplitudes  an  order  of  magnitude 
greater  than  those  necessary  to  induce  Kelvin  Helmholtz  instabilities. 

For  this  reason,  static  wave  breaking  will  be  considered  a "low 
probability"  candidate  for  the  source  of  ocean  fine  structure. 
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Fig.  8 Inertial  Wave  Stability 
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Transverse  steepening  or  resonance  was  the  third  possible 
mechanism  by  which  inertial  waves  could  form  fine  structure.  If  the 
fine  structure  were  created  by  the  transverse  resonance  of  perpendicular 
waves,  then  the  persistence  time  of  the  fine  structure  would  be  of  the 
order  of  l/u  or,  for  inertial  waves,  a few  hours.  The  fine  structure 
data  base  (Table  I)  illustrates  persistence  times  of  the  order  of  days. 
Thus,  to  explain  fine  structure  with  resonance,  we  would  have  to  assume 
that  both  ovservers  (Refs.  2 and  7)  obtained  their  measurements  an 
integral  number  of  wave  periods  apart.  Hayes  et.  al.  obtained  their 
data  at  28°N.  This  corresponds  to  a wave  period  of  24  hrs.  and  would 
require  three  wave  periods  to  be  consistent  with  observation.  Howe  and 
Tait^^^  worked  at  34°N,  or  22  hr.  wave  periods.  Hence,  to  explain 
"persistence  times"  of  33  hrs.  , we  would  have  to  invoke  a resonance 
phenomenon  occuring  at  the  wave  half  period.  While  this  is  all  possible, 
it  does  not  seem  highly  probable.  In  order  to  explain  fine  structure  on 
the  basis  of  resonance  we  must  invoke  large  wave  amplitudes  (£  ^ 1). 
Inertial  waves  of  that  amplitude  are  dynamically  unstable  and  incapabl-i 
of  undergoing  a rever sable  energy  exchange  with  other  waves.  Hence, 
resonance  should  also  be  considered  to  be  a "low  probability"  candidate. 

The  fourth  proposed  mechanism  is  the  self  steepening  of  a plane 
wave  due  to  a longitudinal  velocity  contained  in  the  higher  order  terms 
of  the  dynamic  equations.  Such  a mechanism  requires  that 


U • Grad  1^0, 
w 

where  is  the  fluid  velocity  vector.  For  longitudinal,  or  acoustic  sound 

waves,  U * Grad  I,  is  an  order  £ term.  Since  inertial  waves  are,  to 
w 
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lowest  order,  transverse  waves,  self  steepening  can  be  no  greater  than 
2 

an  e effect.  This  is  consistent  with  the  self  attenuation  of  surface  waves. 

. (18)  . ~ 2 ~2  ~ 

The  dispersion  equation  yields  a (1  + k a ) term  where  a is  the  wave 

^ ^ 2 ^^2  2 

amplitude.  Since  aw  u /a,  k a becomes  e . Hence,  self  steepening 
-2  P 

requires  £ cycles  to  attenuate  the  wave.  The  maximum  value  of  £ that 

we  can  use  to  estimate  the  steepening  time  in  £ . because  all  waves 

c rit 

corresponding  to  £ > £ . are  unstable.  The  steepening  time  T , becomes 

crit  s 


T 

s 


2TT  / 


or 


2TTQ^ 


Estimates  of  the  self  steepening  times  are  illustrated  in  Fig.  9.  Values  of 

T as  large  as  100  days  are  required  in  order  to  be  consistent  with  data, 
s 

It  is  difficult  to  believe  than  an  inertial  wave  would  not  encounter  several 
first  order  effects  during  this  time  period.  These  encounters  would  tend 
to  destroy  any  coherence  of  the  postulated  mechanism.  Hence,  self  steep- 
ening may  also  be  considered  "low  probability". 

Of  the  four  mechanism  proposed,  dynamic  wave  breaking  is  the  most 
probable  mechanism  by  which  inertial  waves  form  fine  structure.  Jhe 
wave  length  scales  and  the  wave  stability  criterion  are  consistent  with 
observation.  However,  it  still  remains  to  be  demonstrated  that  a dynamically 
unstable  inertial  wave  actually  forms  fine  structure. 


IV.  CONCLUSIONS 


The  preceding  analysis  has  made  an  apparent  correlation  between 
ocean  fine  structure  and  inertial  waves.  A priori,  there  are  four  possible 
nonlinear  mechanisms  that  could  form  ocean  fine  structure  from  inertial 
waves:  dynamic  breaking  (shear  instabilities),  static  breaking  (inverted 
density  profile),  transverse  steepening  (resonance)  and  longitudinal 
steepening  (self  steepening  of  waves  riding  on  their  own  velocity  fields). 

Each  mechanism  has  been  assessed  on  the  basis  of  the  available  data.  Of 
the  four  mechanisms  proposed,  dynamic  breaking  is  the  most  probable 
mechanism  by  which  inertial  waves  form  fine  structure.  The  length  scales 
of  the  wave  and  the  stability  criterion  are  consistent  with  observation. 

The  static  breaking  of  internal  waves  would  require  much  larger  wave 
amplitudes  than  does  dynamic  breaking.  Hence,  the  concept  of  static 
instabilities  becomes  academic  since  these  waves  are  already  dynamically 
unstable.  Resonance  phenomenon  would  also  require  larger  wave  amplitudes 
than  does  the  postulate  of  dynamic  breaking.  Thus,  resonance  cannot  occur 
because  inertial  waves  of  the  required  amplitude  are  dynamically  unstable. 
The  possibility  of  self  steepening  has  also  been  assessed.  Estimates  of 
the  self  steepening  times  for  stable  (low  amplitude)  inertial  waves  are  of 
the  order  of  10^  days.  This  is  clearly  a long  time  for  an  inertial  wave  to 
avoid  other  first  order  effects  which  could  destroy  the  coherence  of  the 
mechanism.  Hence,  self  steepening  is  a weak  possibility  compared  to 
that  of  dynamic  wave  breaking. 

Previous  investigators  have  proposed  that  shear  instabilities  lead 
to  ocean  fine  struc<"re.  However,  their  concept  was  that  of  "Billow  Turbu- 
lence" driven  by  gravity  waves  which  propagate  horizontally  along  density 
"sheets"  (discontinunities  in  the  vertical  density  profile).  These  gravity 
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waves  break,  rnix  and  form  the  fine  structure.  The  mechanism  is  illustrated 
in  Figo  10a,  Figure  10b  illustrates  the  present  concept.  A vertically 
propagating  inertial  wave  generates  a horizontal  shear  layer.  The  shear 
layer  is  dynamically  unstable  and  appears  to  break,  forming  the  fine 
structure.  The  wave  length  scales  and  stability  criterion  are  consistent 
with  fine  structure  observations.  However,  it  still  remains  to  be  demon- 
strated that  a dynamically  unstable  inertial  wave  actually  forms  fine 
structure. 
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appendix  a 

INTERNAL  WAVE  MODEL 


Consider  an  orthogonal  coordinate  system  on  the  surface  of  a 
otating  earth  as  shown  in  Fig.  3.  The  rotation  vector  is  denoted  by 
where  ~ denotes  a dimensional  quantity.  The  coordinate  system  is 


chosen  s 


uch  that  n has  no  component  in  the  y direction.  Hence, 


A 


n = i + 0^  k , 


where 


n 


2TT  radians 
12  hours 


(A-1) 


The  low  Mach  number  motion  associated  with  the  ocean  may  be 
treated  as  incompressible.  Hence,  the  Lagrangian  derivative  of  density 

is  zero. 


DP 

D~ 


= 0 


(A- 2) 


The  equation  of  continuity  then  simplifies  to 


where 


\7  • U = 0 , 


^ <-s/  A ^ A 

U=ui  + vj+wk, 


(A-3) 


(A-4) 


The  momentum  equation  includes  the  influence  of  the  Coriolis  forces: 


^ ^ 1 


A 2 


^ + n X u + (u  • V)  u + -I  vp  + gi"  = “u, 

dt  ^ 


(A- 5) 


A-1 


where  V is  the  kinematic  viscosity  and  g is  the  acceleration  of  gravity. 

Equations  (A-2),  (A-S)  and  (A-5)  represent  three  equations  in  the 
density  p,  pressure  p and  the  velocity  U.  To  obtain  a solution,  we  expand  the 
dependent  variables  into  a static  and  wave  contribution: 


P - pQ  + Pg  (z)  + (x,  y.  z,  t).  (A-6) 

P = Po  - Po  g ^ + Ps  + Pw  (A- 7) 

^ ^ ^ /W 

u = (x,  y,  z.  t),  (A_3) 

V = (x.  y,  z.  t),  (A. 9) 

w = w^  (x,  y,  z.  7),  (^_10) 


where  we  have  restricted  the  analysis  to  zero  mean  shear.  The  only  veloci- 
ties are  due  to  the  wave  contribution. 

Setting  all  wave  contributions  equal  to  zero,  we  obtain  the  hydrostatic 

relationship  between  p and  ^ , 

s s 


Pg(z)  g. 


The  density  profile  in  a static  sea  will  be  specified  "a  priori".  Hence,  we 
assert  that 


where  p^  is  the  reference  density  at  z = 0.  The  value  of  p is  expressed 
in  units  of  inverse  length  and  has  values  between  10"^/m  in  the  summer 
thermocline  and  10  /m  in  the  deep  sea. 


A-2 


Substituting  Eqs.  (A-6)  through  (A-10)  into  (A-2),  (A-3)  and  (A-5), 
we  determine  the  equations  of  motion  for  the  wave  contribution; 


dp 

w 

dT 


w 

w 


/Si/  /^ 

r ) p 


w 


= 0. 


(A -11) 


(A -12) 


(A-13) 


where  [i  is  the  viscosity. 

We  define  the  reference  lengths  as  the  inverse  wave  numbers  in  the 
horizontal  and  vertical  directions: 

-1 

^ref  ^ref  2n  H 

and 


The  horizontal  fluid  velocity  is  normalized  to  u^,  an  arbitrary  velocity 
which  will  be  related  to  the  wave  amplitude. 

/s^  ^ 

U . = V r = U 

ref  ret  p 

/s,/ 

The  continuity  equation  suggests  that  the  reference  velocity  in  the  z direction 


/s,/ 

W 


ref 


- u k /k  , 
p H V 


IS 


The  reference  time  is  the  inverse  frequency,  which  will  be  determined  from 
the  wave  equations. 

t - 1/a 
ref 

The  linear  form  of  Eq.  (A-11)  suggests  the  reference  density. 


D . = u p 6 k /a  k 
ref  p H V 


In  the  vertical  momentum  equation,  it  is  not  clear  whether  the  vertical 
pressure  gradient  should  balance  the  buoyancy,  or  the  Coriolis  force. 
Hence,  we  normalize  the  pressure  by  the  buoyant  terms 


ref 


= Q P, 


u a/k 
p H 


but  introduce  the  ordering  parameter  Q which  must  be  determined. 

The  reference  variables  are  used  to  define  the  following  dimension- 
less parameters: 


e = u k /a 
p H 


A = P/k 


V 


and 


6 = k /k 

V K V 


Re  = p a/i-lk 

wo  V 


(Wave  Amplitude), 

(Static  A p/p  Across  Wave), 
(Vertical/Horizontal  Length  Scale), 
(Wave  Reynolds  Number). 


A .4 


The  frequencies  0 and  Q are  normalized  to  N, 


the  Brunt  Vaisala  fre- 


quency, which  is  assumed  to  be  constant. 


~2 

N = g P 

cp  = a / N 

t =Q/N 


The  equations  of  motion  for  the  wave  contribution  are  non 
dimensionalized.  The  statement  of  incompressibility  becomes 


& p — ► 

^ . w + e (U  • 7p  ) = 0. 

Jj  t W W w 


where 


V=  b -^i+6  T— 

V Px  V Sy  dz 


and 

U 

w 

u 

w 

A 

T 4- 

1 + W K 

1 T 

6 ■'  w 

V 

such  that 

-» 

u • 

w 

11 

O (1)  . 

(A-14) 


The  equation  of  continuity  is  expressed  as 


V • U = 0 , 
w 


(A-15) 


and,  expanding  p. 


_P_  = = 1 - A z + » 
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the  X,  y,  and  z momentum  equations  become,  respectively, 


du 


cp 


w 


+ e cp  (U 

ot  w 


V) 


w 


O)o  V 

3 w 


+ cp  Q 


dx 


(A-16) 


cp 


w 


dt 


cp  (U 


V 


w 


w 


0)  u 
3 w 


6 w w 
V 1 w 


+ cp  Q 


V , 

w 


and 


7 ? ? 2 

cp  6 3-^+e6  cp  (U  -r)w  +w  6 cp  V 

^--dt  v w wlv 


V 


w 


2 .2  2 / cp“  \ 2 

+ Q cp  -“+6  p =6  - — r w . 

dz  V w V \Re  / w 


w> 


(A-17) 


(A-18) 


Linear  Solution 

In  the  limit  of  e « 1,  A«  1 and  Re  » 1,  Eqs.  (A-14)  through  (A-18) 

(13) 

are  a linear  system.  From  Phillips  we  know  that  possesses  a wave 
solution  where  the  frequency  is  given  by 

a^=  6 ^/(l  + 6 ^)  + (n-  p)^, 

V V 

A -6 


(A-19) 


where  (-lisa  unit  vector  in  the  direction  of  the  wave  propagation, 

A A A A 

(i=6T6  i + 6 6 j+k, 

1 V Z V 

where 


and 

The  linear  wave  solution  for  w is 

w 

w = A exp  (i  ^ ), 
w 

where  A is  an  arbitrary  constant  and  I is  given  by 

^ = t + y + z* 

From  Eq.  (A-14),  the  density  is  obtained  from  the  vertical  velocity, 

p = - i A exp  (i  I ). 
w 

In  the  limit  of  6 <1,  Eq.  (A-19)  yields  the  nondimensional  frequency, 

2 2 A 7 

cp  =6  + (ou  • (-1)  . (A-20) 

V 

Eq.  (A-20)  indicates  that  there  are  two  limits  in  which  the  momentum  equa- 
tions must  be  solved  for  u , v , and  p . Those  limits  depend  on  the  relative 

w w w 

— A 

size  of  6 and  (uu  • (i),  the  buoyant  and  inertial  frequencies  respectively. 


6 = k /k 
1 X H 


6 = k /k  , 
2 y H 


2 2 
^ =1- 
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For  Gravity  Waves  (6  > (UJ  ' t-L)),  the  solutions  become 


cp  = 6 . 

V 


Q = 1. 


p = A exp  (i  ^), 
w 


u = - 6,  A exp  (i  1), 
w 1 


V = - 6-  A exp  (i^). 
w Z 


For  Inertial  Waves  (((JU  • |a)  > 6^),  the  solutions  are 


Q « 1, 


cp-  0^3. 


u = i (6^  + i 6J  A exp  (i  1), 
w Z 1 


V = - ^0 
w 


(62  + i 6j)  A exp  (i^). 


The  inertial  solution  of  the  vertical  momentum  equation  yields  dif- 

2 

ferent  results  for  p depending  upon  the  relative  size  of  6 and  Cp  . 

w '' 

Physically,  the  different  solutions  arise  due  to  the  relative  size  of  the 
buoyant  and  Coriolis  forces. 

For  Inertial  Gravity  Waves  (buoyancy  balances  pressure),  we  obtain 


Cp  <6  < cp, 

V 

,2,2 

Q = 6 /cp  , 

V 


p -A  exp  (i  ^). 
w 


For  Pure  Inertial  Waves  (Coriolis  force  balances  pressure),  we  obtain  ^ 

2 ; 

6 < cp  , 

V ! 

Q = 6 , “ 1 

V I 

+ i 6i)  A exp  (i  ^). 

The  regimes  of  (p,  b ) space  appropriate  to  gravity  waves,  inertial  gravity  . 

waves  and  pure  inertial  waves  are  illustrated  in  Fig.  (A-1).  The  solutions 
appropriate  to  each  regime  were  given  above. 
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-VERTICAL  WAVELENGTH/HORIZONTAL  WAVELENGTH 


Fig.  A-1  Internal  Wave  Regimes 
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